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Abstract 

The full classical Dirac-Maxwell equations are considered in a somewhat novel form 
and under various simplifying assumptions. A reduction of the equations is performed in 
the case when the Dirac field is static. A further reduction of the equations is made under 
the assumption of spherical symmetry. These static spherically symmetric equations are 
examined in some detail and a numerical solution presented. Some surprising results 
emerge from this investigation: 

• Spherical symmetry necessitates the existence of a magnetic monopole. 

• There exists a uniquely defined solution, determined only by the demand that the 
solution be analytic at infinity. 

• The equations describe highly compact objects with an inner onion like shell struc- 
ture. 



INTRODUCTION 



It is an interesting exercise to compare the current development of a quantum theory 
of gravitation, from the fully non-linear Einstein theory, to the development of QED from 
the linearized Dirac-Maxwell theory. The most startling difference is the large body of 
work on the classical, non-linear, theory of gravitation (General Relativity) - a theory 
which includes, in a self consistent manner, the interactions of the gravitational field 
itself. There is no comparable body of work on the full Dirac-Maxwell theory (Dirac 
equations with electromagnetic interaction. Maxwell equations with Dirac field source - 
the so-called "self interaction"). Of course this situation arose, historically, because of the 
rapid development and stunning success of QED. 

Einstein's equations provide a much better description of gravity than do the linear 
spin-2 equations. Indeed, one can "derive" classical general relativity from the linear, 
massless spin-2 theory by summing all the Feynman diagrams to tree level - see 
and P]. The full Dirac-Maxwell equations should provide a much better description of 
electronic matter than their linearized counterparts (in which self terms are ignored). 

Solutions to the Dirac-Maxwell equations are rare, see |^. There are , however, a 
number of solutions to the Yang-Mills-Dirac and Yang-Mills-Dirac-Higgs equations, see 
1^ and the comprehensive list of references contained therein. 

The paper is organized as follows: in §1 we write the equations in 2-spinor form, this 
description then allows us to (covariantly) solve the Dirac equations for the electromag- 
netic potential and so write down a complete set of equations in terms of the Dirac field. 
In §2 we examine the static and spherically symmetric reductions of the equations. In §3 
we deal with some general properties of the static spherically symmetric system and in 
§4 present a numerical solution to this system. 

§1 The Dirac-Maxwell Equations 

In standard notation the Dirac-Maxwell equations are 




-A, 



(1) 



Q/3 = -4vre 



Employing the 75-diagonal or van der Waerden description, see for example we have 
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with cr^^ the van der Waerden symbols, i.e. 



f^^) = 4 r ° 1 and 



AAl 



V2 I 1 



1 

{a\^ = — X Pauh Matrix, j = 1, 2, 3. 

Where A, i? = 0, 1 and A, 5 = 0, 1 are two-spinor indices (see [Pe]). 
The Dirac bispinor, i\} is 

= \ ^ \ and 

^ = (v^^u^y 

So that the Dirac equations become 

V2 

(2) (5^^ + ,eA'^>^ + ^M^ = 

v2 

where d^^ = = a'^^^A^. 

The Maxwell equations are 

(3) d°'Fai3 = -4:TTejp = -Aire a^^iuAUx + vaVa) 

In the linearized theory the "self current" jp is ignored. We now eliminate the potential 
A^"^ using (2) - another approach is to eliminate "A^*^^^ "using the formal Green's function 
1^ - we will use purely algebraic methods. From the equations (2) we have 

v2 

(4) =ie[A^\v^UB-u\B)] 

However, because of the 2- dimensionality of the 2-spinor space we have 

VaUb - UaVb = eAB{u^vc). 

Here, £oi = = 1, £io = = -1, = £oo = = £n = 0; we define = e^^^b 
and ^A = £ba^^- 

We assume that u'^vc 7^ almost everywhere. Now, = I'W'^fyiP, so u'-^'vc = 
implies that the current vector, j, is null - a massive neutrino field. 
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We can now solve (4) for the electromagnetic potential A, 

(5) A^^ = ^^^LaqBA^^_^^AqbA^^^!!^^^a^A^^a-A^ 

G\U Vc) \^ V 

Notice that, from (5), under the gauge transformation 

— > e*^ I ^^ I we have 

Aa ^ Ace H 9ceX 

e 

as we should expect! 

The four complex equations (2) actually over determine the four real quantities A^. 
We must impose on (5) the condition that A^ is real. These reahty conditions can be 
written as 



(A^^uaUa) = A^^uaUa 

(A^^VaVa) = A^^AVA 



{A^^uaVa) = A^\aua 
With the use of (5) these reality conditions become 

d^^ivAVA) = ^("^^C - ^'Vc) 

(6) UAd^^VA - VaQ^^ua = 

These equations constitute four real first order equations for the four complex quantities 
and v"^. A further four real third order equations for these quantities is obtained 
upon substitution of (5) into the Maxwell equations (3). Note that adding the first two 
equations of (6) leads to the equation of conservation for j". 

§2 Reduction of the System 

2.1 The Static Equations 

Firstly, we impose the condition that the field is static. We assume that there exists 
a Cartesian Lorentz frame in which j"' — S^j^. Imposing this condition one quickly finds 
that 

^ e'^^V^cr^^^M^, with % a real function. 

the current vector is now 
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Now write 

The reality conditions are 

9^i(CV) = ^(ICT + ICT)sinx 
(5o6 + 5ii)(ICT + ICT) = 

C°(5o6 + d,^C - C^(5o6 + 5ii)C° = ^[C°C^(5o6 " ^li) + iCfd,^ " {Cfd,,]x 

the expressions for the potential A^^ can now be written down, although we won't do 
this at this stage. 

Now under a gauge transformation we have 

^ e'l'C^ and 

e 

We fix the gauge by defining real functions X, Y and rj as follows 

^1 — Y e^^^~'^^ 

Our equations can be given in a particularly suggestive three vector form by writing 
(in our Cartesian coordinates) 



with 



and 

y = (2Xycos77, 2Xysin77,X2-y2) 
The reality conditions become 

— + (Vx)xy = o 
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With electromagnetic potential 



o_m _ {X--Y^) dr] (V.v).V 

A - ^ (cos X 1) + 2 e(X2 + y2) dt + 2e(X2 + F^) 



1 



2e(X2 + y2) 



dt 



V + {X^- Y^)Vr] - V X V 



where A^{A\A\A-^) 

The full system is given by the above two sets of equations and the Maxwell equations. 
2.2 Spherical Symmetry 

Wc now impose spherical symmetry upon our static system. A minimal requirement 
that the Dirac field be static and spherically symmetric (in any gauge) is that the vector 
i, above, is spherically symmetric. We require 

[X^,i]=0 

where the Xi, i — 1,2, 3, are the three (vector) generators of rotations. 

These conditions imply that I has time and radial components only and that these 



components are functions of {t,r) only, r = \J{x^Y + (a;2)2 + (a;^)2. 
Using the notation above we have 

+ y2 = R{r) , only 

V = \V\r^Rr 

where f = (sin 9 cos (f), sin 9 sin 0, cos 9) in terms of the polar coordinates r, 9, (f). We have 



X = VRcos{9/2) 
Y = V^sin(^/2) 
77 = (/) 



The Dirac bispinor is now 



i) = e"^"^* ^/R 



The equations are now as follows 



/ _ei(x-^)sin(f) \ 
ei(x+</>) cos(f) 
_e^(x+-^) sin(f) 
\^ e^^^'-^^cosd) j 



x^x{^), R = R{r) 
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1 cot e ^, 

^=7^ (t> 

2 e r 

A =-(cosx-l) + 2^^ 

^{r'^R) = -2mr^i?sinx 
ar 

(7) ^(r^^) = -4V2 7rer2i? 

ar ar 

The really surprising result here is the unavoidable appearance of the magnetic monopole 
term 

1 cot e 

A = (f>. 

Ze r 

Here 4> is the usual azimuthal unit vector, in terms of a coordinate basis A = cos 9 d(f). 

We should also impose a normalisation condition (or finite total charge condition) on 
any solution 

' fdSa < oo 

on any space like hypersurface. This leads to the condition 

r^Rdr < oo 



/■ 



lr=0 

To end this section we rewrite the determining radial equations in a more transparent 
form by introducing the following new (dimensionless) variables 

p — 2mr 
a = 

m 

q = 47rV2—r'^R. 
m 



dy 

— — a + 1 — cos X 
dp 

d ^ ^da 
dp dp 

(8) — = -qsmx 



dp 



§3 Static Spherical Symmetry 

3.1 Some General Properties 



The system of equations (8) possesses the discrete symmetry 



a + 1 ^ -(« + 1) 



This is just the operation of charge conjugation; q needs to be reinterpreted (it was 
originally defined as non-negative) to account for the change in sign of the charge which 
manifests itself on the right hand side of the second equation of (8). We write q — eQ, 
with — 1 and Q > 0, so Q — V^^r^i?; then, under (9), we have e — > — e and Q ^ Q. 
Our equations read, 

dx 

— = a + 1 — cos X 
dp 

dp 

(10) f = -«^'"^ 

Where we have introduced the new variable / to give a set of four first order, ordinary 
differential equations. This new variable is directly related to the magnitude total (Dirac 
field) charge contained in a ball of radius r, 5(r), 



: / fdSa = e / fd^X 
JBir) JB(r) 

= 47reV2 /' s^Ris) ds 

Js=0 



= ^ r Qda 

e Ja=0 

= ^(/(p)-/(o)) 

In view of this - and the fact that ^ — Q — eq is proportional to the charge density on 
a shell of radius r - we will use the following condition on our system (10). 



(CI) 



On p > , / is a bounded function, with bounded first derivative. 

Both / and — have well defined limits as p — > oo 
dp 



We will now develop some qualitative results which indicate the types of solution which 
can exist under rather general (and physically reasonable) conditions. 



Lemma 1 Suppose (%, a, /, Q) is a solution of (10) on p > , then under CI the function 
a is on p > with a and ^ bounded on intervals p > pi > and pa bounded on 

dp ^ 



< p < pi < oo. ///(O) y^O or f (0) = g(0) > then a is unbounded as p — > 0. 



Proof. We first establish that Q and / have well defined limits as p approaches 0. We 
are assuming that the solution (x, a, f, Q) exists on p > 0, so for p2 > pi > and using 
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equations (10), we have 



P2 



|Q(p2)-g(pi)| = I / Q{a)smx{a)da\ 

rp2 

< / \Q sin x\da 
J pi 

< Ml (p2 — Pi) : here Mi — sup Q < oo. 

Letting pi and p2 approach zero, we have \Q{p2) — Q{pi) \ 0. Using Cauchy's criterion 

we conclude that Q{0'^) exists. A similar argument can be given to demonstrate the 
existence of /(O"^). 

The boundedness of ^, on p > pi follows from 

da 2 
% = -^/^ 

Integrate this expression to bound a on p > pi > 0. 
Write il — epa, then 

p— n = — /, and 

dp 



d^Q 
"d^ 

Prom the second of these equations we have, on (0, pi]. 



, , , , Pi , dfl , , dQ , . dQ , , 
So we have (note as a is away from p = so is fi) 

dn 

p— > 0, asp — > 

dp 

Hence from the first of the fl equations we see that fl has a well defined limit as p ^ 0, 
in fact 

limQ(p)=/(0) 

An immediate consequence is that if /(O) 7^ then a is unbounded as p — > 0. 
Now from (10) we have Q > (3(0)e^^ so that /(p) - /(O) > f (0)(1 - e"'') - recall 



dp 

^ = eq = Q. Given the earher result we may take /(O) = - otherwise a is unbounded. We 



can now integrate our ^ equation with this bound for /, 



df 

e(a(p)-a(pi))>^(0) 



(1-e-^) (l-e-^'i; 



Pi 



+ 3^(0) / da 



dp Ja=p (7 



As p — > the integral on the right side of the inequahty diverges to +00. 
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Lemma 2 Suppose {x,a,f,Q) is a solution of (10) on oo > p > 0, then under CI the 
function x is with ^ bounded on intervals p > pi > 0. ///(O) 7^ then x is unbounded 
as p 0. 

Proof. The regularity of x ^ on its presumed interval of existence - is established using 
standard theory (see, for example, [0] or ^) after first noting that the right side of 

dx 

— = a + 1 — cos X 
dp 

is in p (treating known function and using lemma 1) and C°° in x- 

From the above equation we also have 

e(a + l)-l <e^ < e(a + 1) + 1 

which gives the required bounds (using lemma 1). Working on (0,pi) we have (as in the 
proof of lemma 1) 

p dp p 
where Cq = e(a(pi) + 1) - 1 - ^ and = e(a(pi) + !) + !- 
Integrating, 

C2 - C^{pi -p) + fipi) Hp) <ex<Cs- Co(pi - p) + /(O) ln(p) 

with C2 = ex(pi) - /(pi)ln(pi) and C3 = ex(pi) - /(O)ln(pi). If /(O) < choose pi 
near so that /(pi) < and we have from our last inequality that ex ^ 00 as p ^ 0. If 
/(O) > our inequality yields ex — > —00 as p ^ 0. 

There is one other condition which makes sense "physically": if we have an isolated 
system we expect the charge density should go to zero at infinity. 

df 



(C2) 



Q ^ 0, asp ^ 00. 
dp 



Lemma 3 Suppose {x,a,f,Q) is a solution of (10) on p > under conditions CI and 
C2. Then —2 < a^o < where a —>■ as p —>■ 00. 

Proof. We first establish that a has a well defined limit, written as Ooo, as p — 00; for 
P2 > pi > 0, we have 



\a{p2) - a(pi)| 



P2 



fia)da 



< M2(— -—), where M2 = sup |/|. 

Pi P2 R+ 
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Letting pi and p2 approach oo, we conclude that the hmit a — >■ a^x, exists. 
From the first and second equations of (10) we have 

— In a + 1 — COSY = (— + sin V— )/ a + 1 — cosy) 
dp dp dp 

f 

= sm X - — — r 

p^(a + 1 — cos X) 

Consequently, on p > pi > 0, integrating the last equation of (10) we have 

fp 

-ln[Q/Q{pi)] = / sinxia) da 
J pi 

cos Y 

+ 



In 



pi 

a + 1 — cos X 



p f{(y)do 



Pi [a'^{a{a) + 1 - cosx(a-))] 



a(pi) + 1 - cosx(pi) 

Now assume a^o > 0. In fact, this also takes care of the case a^o < —2, since under the 
discrete conjugation symmetry, (9), a ^ — (a + 2) and in particular a^o — > — (ctoo + 2). 
Working on p > pi > we have, for pi large enough, 

a(p) > ao > , for p > pi and some constant ao- 

Then, as a < g < a + 2, so 

< ao < — < M , where Mis a finite constant - see lemma 1. 
dp 

i.e. , ao < a + 1 — cos x < M 

Clearly both terms on the right side of our equation for — In [Q/Q{pi)] are bounded. This 
contradicts our assumption C2 that Q as p — oo. 

The constant a^o can be removed from the potential via a gauge transformation. Under 
■0 —>■ e~*™'*°°*'?/), we have a —>■ a — a^o . After this transformation the Dirac field has 
time dependence e~*^*, where —m < E = (1 + aoo)m < m. 

The three lemmas give a basic characterisation of the solutions obeying CI and C2. 
If /(O) 7^ or Q{0) > then a diverges at the origin, these are solutions which can 
be pictured as an Dirac field surrounding a central charged monopole - the numerical 
solution of §4 is of this type. There is also the possibility of solutions with a and R 
(recall: Q = \pl ^r'^R) everywhere bounded, such solutions were suggested by the work 
of Wakano 0, who examined numerical solutions for what could be called "half linearised" 
Dirac Maxwell equations - "half linearised": if the electrostatic potential is "dominant" 
ignore the Maxwell equation involving the electromagnetic vector potential and vice versa. 
In fact, as the following theorem demonstrates, no such solutions exist. 



10 



Theorem, r/zere does not exist a non-trivial solution of (10) on p > under conditions 
CI and C2 with a and P = Q/p^ bounded on p > 0. 

Proof. From lemma 1 we have, /(O) — Q{0) — 0. Next we estabhsh that both x and a 
have well defined hmits as p — > 0. Note that ^ — a + 1 — cos x is bounded, under the 
hypothesis of the theorem, as p — > 0; hence, by an argument of the sort used previously, 
X ~^ x(0) — Xo, say, as p — > 0. We also have 

f=r Q{a)da= T a'^P{a)da>0 

Ja=0 Ja=0 

Using the mean value theorem we have, for some pi, p > pi > 

f(p) = ppIP{pi) < p^Ms , where Mg = supP < oo 



IR+ 



Thus, 



\a{p2) - a{pi)\ = 



P2 



Pi 



f{p) dp 



< 



rp2 1 

/ pMsdp=-M3{pl-pl) 

J pi z 



Letting pi, P2 — establishes the existence of the limit for a, we write 



lim a — Qq 

p^O 



Now we use an argument similar to that used in the proof of lemma 3 to show that 

dx 
dp 



^ = a + 1 — cos X ^ 0, as p ^ 0. On (0, pi) we have 



[*) - \n[Q{p,)/Q{p)] = In 



a(pi + 1 - cosx(pi) 



a + 1 — cos X 



Ja=p a^[a{a) + 1 



/(cr) da 



cosx((t)] 



Assume, ^(0) = Oq + 1 - cosxo 7^ 0. Then, choosing pi near so that ^(p)^(O) > 0, 



dp ' dp 



for pi > p > 0, we have 



Jo a^aia) + 1 



f{a)da 



[a((7) + l-cosx((j)] 



^1 ada 



dXi 
dp^ 



The right side of this inequality is bounded as p ^ 0. Consequently, the right side of (*) is 
bounded as p ^ 0. But this contradicts the assumption of the theorem that Q — p^P 
as p ^ 0. Thus 



Oo + 1 — cosxo = 0, or 
Oo = -1 + cosxo 
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We now assume e = +1, the case e = — 1 can (of course!) be obtained by conjugation. As 
To = -f/p' < on p > 0, so 



Co = — 1 + cosxo > a > Ooo > — 2 



on p > 0. 

Define new variables 



U = cos(x/2), and 
V = Jq sm{x/2). 



We liave 

dU 1, 

— = — (a + 2)V, and 

dp 2 

dV 1 
—r- = T^o.U. 
dp 2 

The pair U and V also satisfy the following linear, second order equations 

d^U f dU I , ^ 

(n, ^,JL^^i(.,2)v..o 

ap^ p"^a dp 4 

We note that ^;2(^^is bounded on intervals [0,p2], with < p2 < cxd and that ^ is 
bounded on (0, oo] - in the first case we may have Ooo = —2, whereas we may have ao = 
in the second; we also have 

^a(a + 2) < 0, on,p > 0. 

From the definitions of U and V we have f/, \^ — ^ as p ^ or oo. Thus, by the 
maximum principle for odes (see [0)5 conclude U = V = Q,soQ = Q. There do not 
exist non-trivial solutions. 

§4 Numerical Solutions 

Numerical solutions to the system (10), with e = 1, were sought by first expanding in 
a power series from either s = (p = oo, with s = l/p)orp = and then evolving the 



system in s or p, respectively, using a MATLAB interface to the NAG library, |11 
4.1 Solutions Near p = 

In lemma 1 we found that Q = pa was bounded, near p = 0, with 

lim(](p)=/(0) 
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It is natural then to seek solutions of the form a = Q{p)/p, near p = 0, with Q analytic 
in p. From (10) it can be seen that both / and Q must be analytic, with Q{0) = 0. 

However, with Q analytic and Q{0) = 0, the last equation of (10) can only be satisfied 
with Q = - which implies a = Cq + Ci/p, we will refer to such solutions as trivial. The 
behaviour of the system near p = may be quite complex; from the proof of the theorem 
(see equations (11)), with a = Q/p, we see that near p = the second order equations for 
U and V have indicial equation (see 



This implies that U and V have behaviour 

U or K ~ co'i(p) cos( — ^ hip) + ct;2(p) sin( — ^ Inp), near p = 

4.2 Solutions Near p = oo 

Near p = oo we expect a ~ Ooo + ^ + + ' ' Assuming a is analytic in s = 1/p, 
near s = 0, then implies that /, Q and sinx are also analytic in s, in fact 

^ da 
ds 

td'^a 



ds^ 

2 



sin X = 2s + s 



' d^a/ds^' 



^d^a/ds"^ ) 

With the assumption that x is analytic near s = a uniquely defined power series results 
if we demand that the solution be non-trivial (x does have the freedom to add integer 
multiples of 27r). 

The resulting power series has no free parameters, it is uniquely determined. 
The lower order portion of the power series solution is as follows 

^ 21 520 ^ ' 

+ 7 5096 + 
7 2548 ^ ' 

Using the power series to determine initial conditions it was found that the numerical 
results were very stable for a good range of initial values for s (sq = .000001 to sq = .01), 
the results were somewhat unstable for sq < 0.0000001. The results presented in figures 
1 to 4 were obtained by first shooting from near s = towards p = and then using the 
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final values of this run as initial conditions to shoot from near p — towards s = 0, to 
verify the solution. 

In figures 1 to 4 a, f{r) — /(O) (proportional to the "electron" charge interior to 
a ball radius r) and Q are plotted against the radial distance measured in units of the 
Compton wavelength (i.e. against \p — mr). 

Interpretation: The solution represented in figures 1 to 4 can be thought of as a central, 
charged monopole (point source) , surrounded by an oppositely charged Dirac field - near 
oo the electrostatic potential behaves as — ~ and near r = the potential 

behaves as ~ (where a ^ 5.7037 is the magnitude of the slope of the hne in figure 
5, where a is plotted against l/(mr)). 




radial distance, mr radial distance, mr 



Figure 1: The Angular Variable, x- Figure 2: The Potential, a 



At about one half the Compton wavelength^ from the center there is a shielding effect 

and the Coulomb nature of the central charge becomes apparent. At large distances from 
the center the electrostatic charges "cancel" each other. We can calculate the magnitude 
of the total charge due to the Dirac field (see §3.1), with /(O) calculated numerically 
(/(oo) = 0, in this case), 

'L^'"^'^" = ^(/M-/(o)) 

?5i —11.407391 

2e 

This calculation results in a charge of the same magnitude as the central charge. 

The object is highly compact, with a radius of about a (reduced) Compton wavelength 
- see figures 4 and 6. It has an onion like structure consisting of an infinite series of 
spherical shells - the local maxima of Q occur at points, p — p^, where sin(x(p„)) — 0, 
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however from lemma 2 we see that x must diverge as p — > (in the present case) so there 
will be an infinite number of shells. 




12 3 
radial distance, mr 

Figure 3: The "Charge" Interior 
to i?(r),/-/(0). 




O 1 2 3 

radial distance, mr 

Figure 4: The "Charge" on a 
Shell, radius mr, q. 



We note that, in the present = so the Dirac field has time dependence e 

in the gauge for which a — > as p — > 0. The Dirac field has mass 



^ Lfds. - ^(/M-/(o)) 



Tfl 

11.407391 — 
2e2 
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To finish we calculate the energy-momentum density of the system. The symmetric 
energy-momentum tensor is 



Tap 

rpD 



T^, + t:j , with 



rpem 



-a^ {uaUa^P + Va'^a^p) - (wAtiA,a + VaVa^o)] + eA^ajp) 
These expressions are derived from the Lagrangian 

L = \ {uj^d^^UA - UaO^^Uj, - Vj^d^^VA + VaO^^V, 

-^{uAV^ + uj,v^)+ejaA--^^. 



Notice that in the absence of the electromagnetic field, for a Dirac field with time depen- 
dence e~*^*, the energy density is 



In the present case we have 



-00 



rpem 

-'oo 



jo{m + eAo) 



m 



Q 



g2 47]-r2 



:i + a) 



m 



PI 
2 



da" 



1 

V 



dp / 4p2 



These expressions include terms due to the central Coloumb and magnetic monopole 
fields, they lead to an infinite total energy when integrated over IR^ the energy can, 
however, be regularised by removing these singular terms. If we exclude the energy due 
to the interaction between the electromagnetic and Dirac fields, eA(ai(3), then gives 
the mass density mjo as above - however, this procedure is clearly not gauge invariant. 

Finally, it is perhaps worth mentioning that the highly localised "multi- electron fields" 
described here may in fact have applications to objects described in recent experimental 
work |T^ - "geonium" or "kilo-e" objects - consisting of highly localised (point like, from 
the experimental viewpoint) collections of electrons in atomic traps. 
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^ We assume that the constants appearing in the Dirac equation, i.e. e and m have 
there usual meaning - e the square root of the fine structure constant and m the inverse 
of the (reduced) Compton wavelength. 
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